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We calculate the light pseudoscalar decay constants, f K and /}<-, and their ratio using domain 
wall valence quarks and 2+1 flavors of dynamical staggered quarks. Use of the MILC gauge 
configurations allows us to simulate at several sea quark masses and spatial volumes, and with two 
lattice spacings. We study how well our numerical lattice data for light decay constants and meson 
masses is described by next-to-leading order 51/(3) mixed action chiral perturbation theory and 
explain our strategy for the chiral and continuum extrapolation. Combining our result for fx/ fx 
with experimental measurements of pion and kaon leptonic decays allows a model-independent 
determination of iKsl/IV,,^; we find a preliminary value of |V,,. s |/|V„ t /| = 0.2315(45)(7). 
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1. Introduction 



We calculate the masses and decay constants of light pseudoscalar mesons in unquenched lat- 
tice QCD. These quantities serve as benchmarks for lattice calculations because they are relatively 
simple to compute with all sources of systematic errors under control. They also provide a test 
of chiral perturbation theory (#PT), including the more sophisticated versions that include lattice 
discretization effects. Successful calculation of these simple quantities lends confidence to calcu- 
lations of more complicated quantities such as the kaon bag parameter (Bk) and K — > kk matrix 
elements, which also rely on ^PT-guided extrapolations. Furthermore, the ratio of decay constants 
fa/fx is interesting in its own right because it allows for a model-independent determination of the 
ratio of CKM matrix elements (V^I/lV^I [[I], g]. 

Following the approach of the LHP Collaboration, we use HYP-smeared domain-wall valence 
quarks and staggered sea quarks We use gauge configurations generated by the MILC Collab- 
oration with a 2+ 1 flavor improved staggered action because they are publicly available and span a 
wide range of quark masses, lattice spacings, and volumes These lattices allow for good control 
over chiral and continuum extrapolations for light pseudoscalar quantities using staggered valence 
quarks. We show in this work that we have similar control over these extrapolations using domain- 
wall valence quarks. We account for the effects of staggered sea quarks by using the appropriate 
mixed action #PT (MA^PT) expressions [Q]. The #PT formulas for light pseudoscalar quantities 
in the mixed-action theory are more continuum-like than in the purely staggered case, with fewer 
new parameters. The mixed-action approach is especially powerful when considering more com- 
plicated quantities that are not protected from operator mixing under renormalization, such as Bk- 
In the staggered case, mixing with taste-breaking operators under renormalization creates signifi- 
cant complications [Q]. In the mixed-action case, the chiral symmetry of the valence sector makes 
nonperturbative renormalization as straightforward as in dynamical domain- wall simulations [0]. 
The results presented in this work bolster confidence in our future mixed-action calculation of Bk- 

In Ref. [10] we have performed a strong check of the ability of MA^PT to accurately describe 
discretization effects by investigating the isovector scalar correlator. We find that the MA^PT 
prediction for the two-particle intermediate state (bubble) contribution to the scalar correlator is 
in good quantitative agreement with the numerical lattice data, even though there are large dis- 
cretization effects due to staggered sea quarks. Thus we conclude that MA^PT correctly describes 
the dominant unitarity-violating contributions to mixed-action lattice simulations. Fortunately, in 
the case of most weak-matrix elements, MA^PT predicts that non-analytic unitarity-violating er- 
rors should contribute at the sub-percent level on the MILC lattices that we are using. This fact, 
in conjunction with our successful analysis of the scalar correlator, substantiates the claim that 
unitarity-violating effects in mixed-action lattice simulations can be accounted for and removed to 
recover precise continuum values for weak matrix elements. 



2. Lattice calculation and chiral-continuum extrapolation 

We generate data on the MILC ensembles given in Table |T[ The quantities m' s and m' denote 
the values of the simulated staggered sea quark masses, while the unprimed quantities denote the 
physical masses m s and m = (m u + m c i)/2. We compute domain-wall valence quark propagators 
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Table 1: Simulation parameters of MILC staggered gauge configurations used in this work. "Nominal" 
quark masses are shown. 
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with masses from m 4 /10 — m s ; our lightest pion is « 240 MeV. In order to minimize finite- volume 
effects, we restrict the combination m K L > 4. We do not tune the masses of our valence-valence 
pions to any particular values because one cannot recover unitarity to obtain full QCD at nonzero 
lattice spacing in the mixed-action theory. We instead generate many partially quenched data points 
at two lattice spacings, and use MA^PT to extrapolate to the physical quark masses and continuum. 
The one-loop MA^PT expressions for most quantities of interest are the same as for chiral fermions 
except for the appearance of additive shifts to the sea-sea and valence-sea squared meson masses. 
Because the two new splittings are easy to obtain from spectrum calculations [Q, |lC[ ], they do not 
need to be included as free parameters in chiral fits. 

In order to suppress contamination from pions circling the lattice in the time direction, we use 
symmetric and antisymmetric linear combinations of quark propagators with periodic and antiperi- 
odic boundary conditions in our 2-point correlation functions. The approximate chiral symmetry 
of domain- wall quarks allows us to use the pseudoscalar current to obtain f n , and we extract the 
decay constant using the axial Ward identity: 

A WP \/2(m x + m y + 2m res ) 

IP = n — = TE j (2-1) 

VAww mj 

where Aww, ^wp> and m % come from a simultaneous fit of Coulomb gauge-fixed wall-wall and 
wall-point correlators. This channel is statistically cleaner than the axial current, and is protected 
by a non-renormalization theorem. Because chiral symmetry is not exact in our simulations at 
finite Lj = 16, there are small corrections to the axial Ward identity. We account for the dominant 
corrections due to residual chiral symmetry breaking by including m res in Eq. (2A). We estimate 
the subleading 0{m res m q a 2 ) corrections to the axial Ward identity using our results at two lattice 
spacings and largely remove them in the continuum extrapolation. 

In order to convert lattice quantities into physical units we use the MILC Collaboration's de- 
termination of the scale, r\, where r\ is related to the force between static quarks, r\F(r\) = 
1.0 [O, O]. The ratio r\/a can be calculated precisely on each ensemble from the static quark 
potential. We use the mass-independent prescription for r\ described in Ref. [13]. In order to fix 
the absolute lattice scale, one must compute a physical quantity that can be compared directly to 
experiment; we use the T 2S-1S splitting [|14|] and the most recent MILC determination of f K []l! 
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X 2 /d.o.f. = 0.93, CL = 0.58 



5.8 p- 
5.4 - 

£ 7 
+ 4.6 - 

i : 

£^4.2 - 

IN N 

s : 

C 3.8 - 
3.4 - 

3 1 



X 2 /d.o.f. = 0.95, CL = 0.55 



full QCD 

am' /am' =0.01/0.05 
am' /am' = 0.007/0.05 




_L 



_L 



0.2 



0.18 



'() 0.02 0.04 0.06 

ri (m x + m res ) 



0.08 



k 0.16 



0.14 



0.12 




_J_ 



am'/aml = 0.01/1 
am' /am[ = 0.0 124/0.03 1 " 
am'W = 0.0062/0.031 
am'/am f s = 0.0031/0.031 
MILC 
experiment 

J , I i I 



0.02 0.04 0.06 
n (m x + m ies ) 



0.08 



Figure 1: Light pseudoscalar meson mass-squared (left plot) and decay constant (right plot) versus valence 
quark mass. Full QCD curves are obtained using the ratio of the quark normalization factors from tree-level 
%PT to convert staggered quark masses to domain wall masses in the MA^PT formulas. Only degenerate 
(m x = m y ) points are shown. 

The combination of the T mass-splitting and the continuum-extrapolated r\ value at physical quark 
masses leads to the determination Hj )hys = 0.318(7) fm [15]. The use of f K to set the scale yields 
rfy* = 0.3108(15)(+26) f m [|TJ]. This difference between the two scale determinations leads to a 
systematic error in our decay constants labeled "input n" in Table |[ 

We use the SU (3) MA^PT formulas derived in Ref. [Q] to extrapolate our numerical lattice 
data to the continuum and physical quark masses. The choice of SU (3) #PT is appropriate given 
the parameters of our numerical simulations because our light pion masses range from 240-500 
MeV and are not much lighter than the physical kaon, which is integrated out in SU (2) #PT. 
Furthermore, the largest of the taste-splittings on the coarse lattices is not much smaller than the 
kaon mass [a 2 Aj « (460) 2 MeV 2 ], though on the fine lattices it is about a factor of 3 times smaller 
[a 2 Aj « (280) 2 MeV 2 ]. The statistical errors on mp and fp are ~ 0.5% — 2% for most of our data 
points. It is now well-established that NLO #PT does not describe pseudoscalar masses and decay 
constants to percent-level accuracy at the physical kaon mass, nor is it expected to based on power 
counting. Our data set confirms this picture. In order to get good fits (as measured by the correlated 
# 2 /d.o.f.) to even our low-mass data we must include NNLO analytic terms. The two-loop NNLO 
logarithmic corrections to the partially-quenched #PT have been computed in Ref. [|I6|], but we 
have not yet implemented these terms in our fits because of their extremely complicated form. 
These expressions would also have to be modified to account for the staggered sea sector, though, 
given our experience with the one-loop modifications to the mixed action, this is likely a small 
effect. In the region where the NNLO analytic terms that we have added are important, we expect 
the NNLO logarithms to vary slowly enough that their effect is well approximated by the analytic 
terms. Even so, this somewhat ad hoc treatment gives rise to our largest systematic error in the 
determination of decay constants. 

We perform chiral fits using two different mass ranges: one with all pions less than 500 MeV in 
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X 2 /d.o.f. = 1.20, CL = 0.15 




n (m x + m KS ) r i (m x + m res ) 



Figure 2: High-mass fit to fp data with chiral continuum extrapolation curves for and f n on left. For 
nondegenerate mass points, the lighter valence quark mass is shown along the x-axis. Again, full QCD curves 
are obtained using the ratio of the quark normalization factors from tree-level ^PT to convert staggered quark 
masses to domain wall masses in the MA^PT formulas. Study of convergence of ^PT for decay constants 
on right. 



order to study the chiral expansion, and another including masses up to the strange quark, which is 
used to determine our central values for masses and decay constants. We correct all data points for 
finite-volume effects using one-loop MA^PT, and take as a conservative estimate for the residual 
finite- volume errors the entire one-loop correction at the data point where that correction is largest. 
The light-mass fits to mp/2(m x + m res ) and to fp are shown in Figure []]. For the leading-order low- 
energy constants (LEC's) appearing in the chiral Lagrangian we use the physical values (J n , jj.), 
rather than the values obtained in the SU (3) chiral limit. To one-loop order this is consistent, and 
studies by both the MILC and JLQCD Collaborations suggest that use of a physical parameter for 
the chiral coupling (f K rather than fo) is likely to result in a more convergent chiral expansion [0, 



17]. We fit the pseudoscalar masses and decay constants separately, since they have no common 



parameters through one-loop once the leading order LEC's are set to the physical values. We vary 
the parameter / appearing in the one-loop expression between the values /o and fa as part of our 
estimate of the chiral fit systematic error. Each of our low-mass fits has 42 data points and 9-10 fit 
parameters; our high-mass fits have 69 data points and 15-16 fit parameters. We do not constrain the 
size of the LEC's with priors. Figure [I] shows the continuum full QCD curve for f n with statistical 
errors as a cyan band. The black curve on the decay constant plot is the continuum full QCD result 
determined by MILC using staggered fermions, after rescaling the bare quark mass in the jc-axis 
by the appropriate factor. Despite the very different shapes of the partially quenched mixed-action 
data and staggered data (not shown; see Ref. [f2|]) at finite lattice spacing, the continuum curves are 
in good agreement. We also show the experimentally-measured value of f % (which is offset for 
clarity) using the T spectrum to set the scale for comparison. 

We present the result of the high-mass fit used to obtain f n and fx in Figure ||. The cyan and 
grey bands show the continuum full QCD extrapolations with statistical errors for degenerate (fx) 
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and nondegenerate (fx) decay constants, respectively. For the nondegenerate curve, the strange 
mass is fixed at its physical value. We determine the physical bare valence u-, d-, and s-quark 
masses from our fits to mp; the physical bare sea quark masses were determined by MILC in a 
similar manner. In order to obtain the charged pion and kaon decay constants, we extrapolate the 
light valence-quark mass to m for f„ but to m u for fx- The light sea-quark masses are set to m for 
both decay constants, and this leads to a slight vertical offset in the final value for fx as compared 
to the grey full QCD band in Figure ||. In order to achieve good fits including nondegenerate 
pseudoscalar masses up to « 600 MeV, we must include higher-order polynomial terms. Figure || 
shows all of the data points used in the high-mass fp fit. In addition to the NNLO analytic terms, 
we include terms proportional to (m x + m y ) n . We find that we cannot get good correlated fits unless 
n is as high as six. Figure |2] compares our results for decay constants with those obtained by the 
MILC Collaboration. This comparison is more meaningful than a comparison with experiment 
because we are using the MILC r\ (obtained from f K ) to set the scale, and one needs |V„ S | in order 
to extract fx from experiment. Figure || also shows the breakdown of our high-mass fit including 
terms through LO, NLO, and "all higher orders" to the degenerate SU(3) curve in the continuum 
limit. The right-most part of this plot corresponds to 500 MeV, where we do not expect ^PT to 
be especially convergent. Because we are interpolating in the quark mass in this region, we expect 
the high-order polynomials to approximate higher-order nonanalytic effects, whereas closer to the 
physical pion mass, the NLO contributions should be dominant. We find the results of our fits to be 
consistent with the expectations from chiral power-counting. It would, nevertheless, be valuable to 
continue this study with the complete and correct NNLO formula including chiral logarithms. 

The error in Table ^| labeled "chiral-continuum extrapolation" is estimated by performing a 
number of different "reasonable" fits and taking the spread between them. For example, we add 
additional higher-order analytic terms (giving fits with acceptable confidence levels) and compare 
the results with those obtained from the preferred fit. We also vary the coefficient of the log terms 
between /o and f n as part of our estimate of the f„ error, and between f n and fx as part of our 
estimate of the error in fx- Although we include terms proportional to a 2 in the preferred fit, 
there is some ambiguity (with only two lattice spacings) in the dominant source of discretization 
errors, which may be purely a corrections, taste-breaking terms proportional to a 2 a 2 , or chiral 
symmetry breaking terms proportional to m res a . It is also possible that the different sources lead 
to discretization effects of the same size. We therefore vary the coefficient of the NLO a 2 analytic 
terms and include the resulting spread in the decay constants as part of the systematic error. We 
also include the parametric uncertainty coming from the uncertainty in the bare quark masses in 
the chiral extrapolation error. 



Table 2: Preliminary error budget. Uncertainties are shown as percentages. 



source 


h 


fn 


fKlf% 


statistics 


1.1 


1.5 


1.3 


input r\ 


1.6 


2.0 


0.3 


chiral-continuum extrapolation 


2.3 


2.2 


1.0 


finite volume 


0.3 


0.9 


0.9 


total error 


3.0 


3.4 


1.9 
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3. Results and Conclusions 

Our preliminary results for the light pseudoscalar meson decay constants and their ratio are 

A = 129.1(19)(40) MeV, f K = 153.9(17)(44) MeV, f K /f n = 1.191(16)(17), 

where the first error is statistical and the second is the sum of systematic errors added in quadrature. 
A breakdown of the total error for each quantity is given in Table ||[ Our errors are somewhat 
reduced as compared to the results presented in the talk, due to the completion of an additional 
fine ensemble with a light strange sea quark. Given the experimental values for the kaon and 



pion leptonic branching fractions and the electroweak corrections [18], we obtain = 
0.2315(45)(7), where the first error is the lattice error (with statistical and systematic errors added 
in quadrature), and the second is the combined error from experiment and electroweak radiative 



corrections. Taking the PDG value of \V ui j\ = 0.97418(27) from superallowed j3 -decay [18], we 



obtain \V US \ = 0.2255(44) (7). This is consistent with the value of \V US \ = 0.2255(19) coming from 



semileptonic kaon decays and non-lattice theory [18]. Given the above value of \V u d\, our value of 
\Vus\ is also consistent with unitarity. 
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